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ABSTRACT 

We present ~ 3" resolution maps of CO, its isotopologues, and HCN from in the 
center of Maffei 2. The J=l-0 rotational lines of ^^CO, ^^CO, C^^O and HCN, and the 
J=2-l lines of ^^CO and C^^O were observed with the OVRO and BIMA arrays. The 
lower opacity CO isotopologues give more reliable constraints on H2 column densities 
and physical conditions than optically thick ^^CO. The J=2-l/l-0 line ratios of the 
isotopologues constrain the bulk of the molecular gas to originate in low excitation, 
subthermal gas. Prom LVG modeling, we infer that the central GMCs have nn^ ^ 
10^'^^ cmr^ and T^ ^ 30 K. Continuum emission at 3.4 mm, 2.7 mm and 1.4 mm was 
mapped to determine the distribution and amount of H II regions and dust. Column 
densities derived from C^'^^O and 1.4 mm dust continuum fluxes indicate the standard 
Galactic conversion factor overestimates the amount of molecular gas in the center of 
Maffei 2 by factors of ~2-4. Gas morphology and the clear "parallelogram" in the 
Position- Velocity diagram shows that molecular gas orbits within the potential of a 
nuclear (^^220 pc) bar. The nuclear bar is distinct from the bar that governs the large 
scale morphology of Maffei 2. Giant molecular clouds in the nucleus are nonspherical and 
have large lincwidths, due to tidal effects. Dense gas and star formation arc concentrated 
at the sites of the xi — X2 orbit intersections of the nuclear bar, suggesting that the 
starburst is dynamically triggered. 
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Introduction 



Concentrations of molecular gas are common in the centers of large spiral galaxies, often in 
the form of nuclear bars. Bars represent likely mechanisms by which rapid angular momentum 



loss and gas inflow concentrate gas in the centers of galaxies (eg. Sakamoto et al. 1999: Sheth et 



al. 2005 Knapen 2005). Secondary nuclear bars can be important in controlling the dynamics of 



the innermost few hundred parsecs of galaxies (eg., Shlosman, Frank Sz Begelman 1989 Friedli Sz 



Martinet |[T993l [Heller et ai:][200]l [Maciejewski et al.||2002| [Shlosman Heller|[2002l [Englmaier El 



Shlosman 2004 1 . In addition to being an avenue by which nuclear gas and stellar masses are built 



up, these bars can also significantly influence the physical and chemical properties of molecular 



clouds within them (eg. Meier Turner 2001 Petitpas Wilson 2003 Meier fc Turner 2004 



20051. 



We have observed the two lowest rotational lines of ^^CO and C^^O, plus new high resolution 
images of ^^CO(l-O) and HCN(l-O) in the nearby galaxy, Maffei 2, with the Owens Valley Millimeter 
(OVRO) Array and the Berkeley-IUinois-Maryland Association Array (BIMA). Because they are 
optically thin, or nearly so, CO isotopologues more directly trace the entire molecular column 
density than does optically thick ^^C^^O (hereafter "CO"). Also ratios amongst isotopologues are 
more sensitive to changes in density and temperature throughout the clouds. HCN, which has a 
higher critical density than CO, constrains the dense molecular cloud component. Here we aim to 
characterize the properties of molecular clouds in the center of this barred galaxy, and to connect 
those properties with the dynamics of the nucleus. 

Maffei 2 is one of the closest large spirals (D~3.3 Mpc, §2; Table [T|, but lies hidden behind 



more than 5 magnitudes of Galactic visual extinction (Maffei 19681. A disturbed, strongly barred 



spiral galaxy (Hurt et al. 1993a; Buta Sz McCall 19991, Maffei 2 has an asymmetric HI disk and 



tidal arms that suggest a recent interaction with a small satellite (Hurt, Turner Sz Ho 19961 



The interaction may be responsible for the bright nuclear CO emission (eg., Rickard, Turner &z 



Palmer 1977; Ishiguro et al 



1989 



Lob — 1.7 X 10 (Turner Sz Ho 1994, corrected for distance). 



Mason Sz Wilson 2004) and active nuclear star formation of 



Distance to Maffei 2 



Galactic extinction complicates distance determinations for this nearby (Vl5_r=-30 km s ^) 
galaxy. Some have suggested that Maffei 2 is close enough to have a significant dynamical influence 



on the Local Group 
the members of the 


(Buta k McCall[[1983 


Zheng, Valtonen & Byrd[[1991 


. Estimated distances to 


IC 342/Maffei 2 group range from 1.7-5.3 Mpc ( 


Buta & McCall[[1983 


McCall 


1989[ Luppino Hz Tonry[[1993t [Karachentsev Hz Tikhonov[[1993 1994[ Krismer, Tully Hz Gioia[[1995l 


Karachentsev et al. 


1997 


Ivanov et al. 


1999 


Davidge Sz van den Bergh 


2001 


). Closer distances 



(Dmpc ~ 2) tend to come from the Faber-Jackson relationship and the brightest supergiants method 
( [Buta k Mc"CaIIl[T983| [Karachentsev Tikhonov||T993l [19941 [Karachentsev et al.[|T997| , while the 
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farther distances (Dmpc ~ 4-5) come from Tully-Fisher relations on Maffei 2 and surface brightness 
fluctuations methods towards its companion Maffei 1 ( |Hurt||1993| [Luppino Tonry||1993| [Krismer, 



TuUy Sz Gioia 1995 Ivanov et al. 1999). In several cases, the same method gives wide ranges of 



distances for different galaxies within the same group, suggesting the group is spatially extended 



(likely) or each distance measurement has higher uncertainties than claimed (Krismer, TuUy &: 



Gioiaj 1995| Karachentsev et al. 


1997 


). Recent studies appear to be convergin 


g to D ~ 3-3.5 Mpc 


for the IC 342/Maffei 2 group ( 


Saha, Claver, Sz Hoessel 


2002 


Fingerhut et al. 


2003 


Karachentsev 


et al. 


2003; Karachentsev 


2005 


Fingerhut et al.||2007 


1. 


Fingerhut et al. 


(2007 


do a self consistent 



analysis of several different measurements and obtain a distance of 3.3 Mpc. We adopt this distance 
for Maffei 2, with uncertainties of ~ 50 %. Quoted uncertainties in this paper do not include this 
systematic uncertainty. 



3. Observations 

Aperture synthesis observations of the 13CO(1-0), 13C0(2-1), 01^0(1-0), 01^0(2-1) and HCN(1- 
0) lines were obtained with the Owens Valley Radio Observatory (OVRO) Millimeter Array between 
1993 October 26 and 1999 March 29. The ^^00(1-0) and ^300(2-1) data were obtained when the 
OVRO array had five 10.4 meter antennas, while the remaining data are from the six-element array 
(Padin et al. 1991 Scoville et al. 1994 1 . Observing parameters are listed in Table [2] Separate 1 



GHz bandwidth continuum channels at 3.4 mm, 2.7 mm and 1.4 mm were also recorded. 3C84 
and 0224-1-671 were used to calibrate instrumental amplitudes and phases. Absolute fluxes were 
calibrated using Neptune, Uranus and 3C273 as standards, with additional observations of 3C84 
and 3C454.3 as consistency checks. Absolute flux calibration should be good to 10 - 15% for the 3 
mm data and 20 - 25% for the 1 mm data, and internally consistent between each transition. 

We have also obtained high resolution (~ 3 ) observations of the ^^CO(l-O) transition with 
the ten element Berkeley-IUinois-Maryland Association (BIMA) Arrayp] ( [Welch et al.||l996 |. Phase 



calibration was done with 0224-1-671 and the ultracompact H II region W3(0H) was used for flux 
calibration. 

Each OVRO track includes at least two phase centers separated by less than the FWHM power 
points of the primary beams of the dishes (Table [2|. The pointings were naturally weighted and 
mosaiced using the MIRIAD. Quoted noise levels are the rms from line-free channels of the spectral 
cube half-way between the map centers and edges. The true noise level is slightly lower (~ 10%) 
at the phase centers and somewhat higher towards the edge of the primary beams due to the 
mosaicing. Subsequent data analysis was done with the NRAO AIPS. 

Since interferometers act as spatial filters, it is important that emission maps being compared 



Operated by the University of California, Berkeley, the University of Illinois and the University of Maryland with 
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have similar (n, v) coverage, including minimum baselines. For the 3 mm lines observed with 
OVRO, the {u, v) coverage is very similar. For the J=2~l lines at A = 1 mm, the (n, v) coverage 
is consistent (C^^O(l-O) and C^*0(2-l) were observed simultaneously), but scaled up by a factor 
of two from their 3 mm counterparts ( {u,v) = Bproj/X, projected baselines in the east-west and 
north-south directions respectively). The 2-1 transitions were tapered to match the {u,v) range 
of the corresponding 1-0 transition. For the OVRO datasets, the minimum (1-0) [(2-1)] baselines 
are {u,v)min — 5.5 k\ [11 A;A]. Thus the images are insensitive to emission on spatial scales ^35- 
40" (~500-600 pc) for 110 GHz and >20" (~300 pc) for 220 GHz. For the ^^00(1-0) transition, 
observed at BIMA, uv coverage is similar to the OVRO 3 mm uv coverage ((u, v)min ~ 5.4 kX). 

We have estimated the amount of flux resolved out of each map due to missing short spacings. 
Single-dish spectra of 1^00(1-0), 13CO(1-0), ^3C0(2-1) and HCN(l-O) exist for Maffei 2 in the 
literature (Table [2j) The interferometer map for each line was convolved to the beamsize of the 
single-dish and then sampled at the same pointing center. Within uncertainties, all of the flux of 
the CO(l-O) and 1^00(1-0) lines is detected. The 1^00(2-1) and HCN(l-O) maps recover ~ 60% 
of the single dish flux, although the single-dish ^^C0(2-l) flux measurement is rather uncertain 
(Wild et al. 1992 1. We expect the interferometer to recover fractions of C^^O flux similar to the 



corresponding ^^CO lines. 

To generate integrated intensity maps, a mask was made by convolving the channel maps to 
lO" resolution, then blanking regions of emission < 2a. This mask was then used to blank out 
non-signal portions for the full resolution channel cube. Velocities from -160 km s~^ to 100 km 
s~^ were integrated, including emission >1.2a in the full resolution channel maps. For the line 
ratio maps, the integrated intensity maps were convolved with an elliptical Gaussian to the beam 
size of the lowest resolution maps (^^CO(l-O) and C^^O(l-O)). This gives a resolution of ~ 3.9" 
for the line ratio maps. Regions of emission < 3cr in either line {2a for the C^^O(2-1)/C^^O(1-0) 
map) were blanked in making the ratio maps. Because of the uncertainties in fluxes and absolute 
positions, the ratio maps are estimated to be accurate to ~ 35-40% in magnitude, and ~ 2" in 
position (excluding possible systematic errors associated with differences in resolved-out flux). 

Since Maffei 2 is within 1° of the Galactic Plane, and at essentially zero redshift, we have 
to consider the possibility of contamination from Galactic CO. Galactic HI emission significantly 



affects VLA images of Maffei 2 in the velocity range -70 km s ^-0 km s ^ (Hurt, Turner Sz Ho 
1996), with the strongest absorption at -40 km s^^. Galactic HI emission is more widespread than 



CO emission, however, and what CO emission there is near Maffei 2 has very narrow lines (< few 
km s~^, based on our spectra from the late NRAO 12 Meter Telescope.) Inspection of our channel 
maps reveals no obvious evidence of Galactic CO emission, and so we disregard it. 
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Molecular Gas in the Center of Maffei 2: Overview 
4.1. Gas Morphology 



CO emission in MafFei 2 takes the form of two prominent and highly inclined arms (Ishiguroj 
et al.|[T989 |, which form the molecular bar, as shown in the integrated intensity maps of Figure [Tj 
The emission extends roughly l!5 (1440 kpc) along the major axis (Table [3|. The brightest CO 
emission emerges from clouds within the central 15" (~240 pc) of the galaxy. 

A higher (2"; ~30 pc) resolution, uniformly-weighted image of CO(l-O) (Figure|2^) shows that 
the bright CO peaks resolve into GMCs. The central two CO peaks become a nuclear ring of radius 
~5", or ~80 pc about the dynamical center. The eastern side of this ring is brighter in CO(l-O) 
but ^^CO(l-O) remains rather uniform. The molecular arms are roughly linear features running 
northeast and southwest, terminating at the central ring. Peak observed brightness temperatures, 
Tmb: reach ~ 31 K, and are typically ^ 10 K across much of the arms. 

Cloud properties — size, linewidth, temperature, mass — are derived for the brightest molecular 
clouds in Maffei 2, using the uniformly-weighted CO(l-O) image (Figure [2|i). Following Meier Sz 



Turner (2001 2004 1 (Table |4|, a molecular cloud is defined as a region of spatially and spectrally 



localized emission greater than 2a in two adjacent channels, but need not necessarily be a grav- 
itationally bound entity. Each cloud was fit from channel maps that include only the gas over 
its localized velocity range. Clumps separated by one beamwidth or less are considered the same 
CMC. Cloud complexes are labeled A-H based on their locations in the lower resolution maps; 
sub-clumps resolved in the higher resolution image are numbered. 

Most of the GMCs are resolved along at least one axis, with sizes of ~40-110 pc (Table |4]). 
The clouds are significantly elongated, with axial ratios often greater than 2. Typical position 
angles of the clouds, 20-60°, are very similar to the sky plane position angle of the bar, ~40°. This 
is not an artifact of the beam shape, since the beam is elongated perpendicular to the bar. Nor 
is the elongation due to an underlying smooth gas component along the bar, since fits come only 
from maps localized in velocity. If the elongation of the clouds is a foreshortening effect due to the 
high inclination of Maffei 2 (67°; Table [T|, then the GMCs must be flattened perpendicular to the 
plane of the galaxy, that is, disk-like, rather than spherical. However, since similar elongations are 



observed in the molecular clouds along the bar in the nucleus of the face-on galaxy IC 342 ([Meier 



k Turner|[200T| ), we consider cloud elongation along the bar more likely. The shapes of the nuclear 



GMCs are clearly affected by their location within the bar. 

With the exception of two small GMCs (D2 and HI), cloud linewidths are > 50 km s~^ FWHM, 
and approach ~ 100 km s~^ in a couple of locations. If these clouds were in virial equilibrium, then 
their individual masses would be in excess of 1O''M0. However, these clouds are very unlikely to 
be in virial equilibrium (^. This is further demonstrated by the fact that there is no correlation 
between the size (Vab) and Avi/2 in Table |4| 
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The CO isotopologues generally follow the brighter CO emission, but there are subtle differ- 
ences. Weak ^^CO(l-O) emission (Fig. Wp) extends to the map's edge, roughly along the major 



axis of the large-scale near-infrared (NIR) (Hurt et al. 1993a) and the large-scale molecular bar 



(Mason & Wilson 2004|. The ^^C0(2-l) emission (Fig. [ip) is found mostly at cloud peaks due 
to the higher critical density {Aij) of the J=2~l line, but some diffuse gas not associated with the 
clumps is resolved out. Little emission off the trailing southwestern CO arm ridges (CMC H) is 
seen in ^^CO. 

C^^O emission (Figs. ^ andjTjf) follows ^^CO, but the C^^O linewidths are slightly narrower. 
This may be a critical density effect, with the lower opacity C^^O more confined to the dense 
cores. As a result of these spatial differences, comparisons of C^'^0 line intensities with ^^CO and 
^■^CO will slightly overestimate their true temperature ratios. Peak main-beam temperatures are 
T^b ~ 0.5 K (T„fc ~ 1.0 K) in a 3'.'9 x 3'.'4 beam (2'.'4 x 2'.'3) for 01^0(1-0) (01*0(2-1)). 

The HCN(l-O) map is shown in Figure [T|i. Because of its much larger electric dipole moment 
(~30 times CO), HCN has a critical density nearly 1000 times higher than CO, and traces high 
density gas. HCN predominately traces the two inner peaks of GMCs D-l-E and F. HCN(l-O) 
emission falls off with distance from the center of the galaxy faster than seen in any of the other 
lines (note particularly CMC G). There is also evidence that the HCN is more strongly confined 
to the GMCs than CO or ^'^CO. Apparently the densest molecular gas is localized more strongly 
to the very center of the galaxy. 



4.2. Star Formation in Maffei 2: Millimeter Continuum Images of H II Regions and 

Dust 



Continuum maps of Maffei 2 at 3.4 mm, 2.7 mm and 1.4 mm are presented in Figure |3] along 



with the 2 cm VLA maps from Turner Sz Ho ( 1994). The 3.4 mm continuum has been corrected for 



the contribution of HCN and HCO"*" lines within the bandpass, and the 2.7 mm continuum map for 
the contribution of ^'^CO. The advantage of imaging continuum at 3 mm is that this is the part of 
the spectrum where the free-free emission component from H II regions is at its maximum relative 
to other sources of emission, such as nonthermal synchrotron and thermal dust emission. 

There are three main 3.4 mm continuum sources near the center of Maffei 2, with weaker 
sources towards the southwestern bar end (GMC G) (Table [5]). Four central sources are found at 2 



cm; source HI has a non-thermal spectral index between 6 cm and 2 cm (Turner Sz Ho 1994), and 



is predictably absent in the millimeter maps. Sources I & II are coincident with GMCs D and E 
and each have fluxes of ~ 5.8 mJy. Source IV is just north of GMC F and has a flux of ~ 5.0 mJy. 
The non-thermal source HI is not associated with any of the bright GMCs. At higher resolution 



these continuum sources resolve into a collection of SNR and H II regions (Figure 2d; Tsai et al. 



2006). 



Spectral energy distributions (SEDs) for each of the main radio sources are shown in Figure |4j 
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Three components are fit: synchrotron, bremsstrahlung (free-free) and dust, with spectral indices 
of -0.7, -0.1 and 3.5 (/3 = 1.5) respectively: 



"^^^44.885 J +"^^^88.92 J I219.3J ' 

Estimated fluxes for free-free and dust emission are recorded in Table [6| At cm wavelengths, the 
central continuum sources have spectral indices between 6 cm and 2 cm, of -0.63 [Sy oc i^"), and 
therefore are dominated by synchrotron emission. The spectral index between 2 cm and 3.4 mm, 
a| 4, flattens to -0.1 - -0.3, as expected for H II regions dominated by free-free emission. There can 
be mixtures of synchrotron and free-free emission within our beam; we have also shown cm-wave 



fluxes for the compact sources (Tsai et al. 2006 corrected for distance) for comparison in Figure 



[4] Our fits indicate that towards I, II and IV the 3 mm continuum emission is dominated by the 
compact, free- free emission sources. 

The 1.4 mm continuum map is shown in Figure [3]i, convolved to the resolution of the 2.7 mm 
map. Emission peaks towards Source II at ~21 mJy beam^^. Continuum fluxes at 1.4 mm are 
larger than the 3.4 or 2.7 mm fluxes, indicating a rising spectral index between 2.7 mm and 1.4 
mm, 02 7' +1-5. The 1.4 mm emission is therefore a mixture of free-free and dust emission, with 
the predominance of dust varying with position. The total flux associated with dust emission, after 
removing the estimated thermal free-free contribution, is S\Amm ~7-19 mJy for each source. 



5. Gas Excitation and Opacity Across the Nucleus of MafFei 2 

Excitation temperatures are important for understanding molecular gas properties and how 
they vary across the nucleus. The J=2-l and J=l-0 lines are sensitive to relatively cool gas in 
GMCs, and the low J CO lines, especially CO(l-O), are thermalized in all but the lowest density 
molecular clouds. CO isotopologues thermalize at somewhat higher densities (^10^ cm~'^) due to 
their lower opacity, making them excellent probes of gas excitation in this density regime. 



5.1. Excitation Temperatures 

Excitation temperatures, Tgx, are constrained by the ratios of integrated intensities of the J 
= 2-1 and J = 1-0 lines, 

f Tiodv ^ ^/lo(Vio(re,) Jl0(Te„^f.))(l " e'^-io) ' 

where ^Ju{Tex) = {hu / k) / {exp{hu / kT^x} — I) , ^Tui and are the optical depth and filling factors 
of the iih isotopologue, respectively. We assume LTE (constant Tex with J) throughout the cloud. 
Limitations of this assumption are noted below. 
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The ^3CO(2-l)/^3cQ(i_o) line ratios for the nuclear bar are shown in Figure H (Table |7|. 
Values range from 0.3 to 0.8, corresponding to Tea; — 3-6 K if ^^CO is optically thin, or up to 10 
K, if completely thick. This ratio peaks towards the central two GMC complexes (GMC D+E &: 
F), and falls with increasing radial distance. C^®O(2-l)/C^^O(l-0) is also higher at the central 
two emission peaks, with values of ~ 0.68-0.79. C^^O is almost certainly optically thin. From 
CiSO(2-l)/CiSO(l-0), we obtain T^x ~5.4-6 K (Tableg. There are several regions off the GMCs 
(between C and D; H2) that have higher ratios. The ratios are largest between clouds: perhaps 
the intercloud gas is warmer than the GMCs (although emission is weak here) . Figure [6] shows the 
average peak T^b ratios schematically as a function of velocity to show changes in excitation along 
the central ring. The eastern side of the central ring has the highest T^x and this is the side closest 
to the starburst. 

Tex based on the isotopologues are significantly lower than both the Tex implied by the single 



dish CO(2-1)/CO(1-0) line ratio of ^2 ( Sargent et al.|[l985 ) and the brightness temperatures, T^b 



estimated from the high resolution CO(l-O). Single-dish CO(3-2)/CO(l-0) line ratios are also 

as are i3c0(3-2)/i3C0(2-l) ratios (1.6; 



high, ~1.3-1.8 ( 


Hurt et al. 


1993b 


Wall et al. 


1993 


. These CO ratios 



Other molecules indicate a range of gas temperatures. From ammonia, Takano et al. (2000) find a 
rotational temperature of Trot ~30 K that is constant across the field, and an ortho-to-para ratio 
consistent with formation at 13 K. Henkel et al. (20001 obtain Trot= 85 K, based predominantly on 
the inclusion of the the high energy metastable transition, (J,K) = (4,4). However, they point out 
that it is possible to fit the four lowest metastable transitions with a cool component and a warm 



component. Rieu et al. (1991 1 derived a low Tg^ of 10 K from a multi-line study of HNCO. 



Excitation of molecular clouds in the nucleus of Maffei 2 is complex, and different molecular 
transitions will find different values for Tex, depending on where the molecules are found. Some of 
the differences in line ratios between CO, its optically thinner isotopologues, and other molecular 
tracers, are due to the isotopologues being subthermally excited relative to CO, such that Tex < T^. 
because < ricritiCO). If the densities determined from the LVG analysis are correct then the 
Tex of the isotopologues imply kinetic temperatures of Tj. = 15 - 35 K (^. These are close to but 
still slightly cooler than the (cool component of) ammonia. Thus the bulk of molecular gas, traced 
by the optically thinner species, is cool. However CO emission is unlikely to be subthermal at these 
densities. So the high single-dish CO line ratios are inconsistent with the physical conditions of 
this component and require the presence of some warmer gas. Whether the emission from the high 



opacity CO transitions originates in warmer envelopes of the clouds (as in IC 342, Turner et al 



1993 Meier &; Turner 2001), or from a compact, dense component (possibly associated with the 
high temperature ammonia component) remains unclear from the current data. The low Tex of 
the higher critical density HNCO may argue against the latter, but spatially-dependent chemical 
effects may also be involved here. 
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5.2. Opacity of the ^^CO and C^^O Lines 



The opacity of the ^'^CO hne is important for mass determinations and the interpretation 
of brightness temperatures. Based on the large CO(3-2)/^^CO(l-0) ratio, Hurt et al. (1993b I 
estimated that the ti2co(i-o) ~ O-^i or that tco{i-o) is only 5-7, on the assumption that CO(3-2) 
and ^3CO(1-0) have the same Tg^. With higher spatial resolution this now appears not to be the 
case. 

Better opacity estimates are obtained by avoiding ratios taken between lines with widely differ- 
ent opacities, particularly in situations where temperature gradients and other non-LTE effects may 
be present. The ^^CO(1-0)/^^CO(1-0) integrated intensity ratio map of Maffei 2 is shown in Figure 
[5^. Values range from 2.4-4.3. The line ratio is lower than expected if the ^^CO(l-O) line (and the 
^^CO(l-O) line) have negligible opacities for adopted abundance ratios of [H2]/[^^CO]=4.7 x 10^ and 
[H2]/[Ci80]=2.9 X 106, or [i^COJ/pa] = 8.5 x 10-^ [i^COJ/pCO] = 60 and [Ci6O]/[Ci8O]=250 
dFrerking et~aI||T982| [Henkel et al][]i94| [Wilson Rood||T994t |Wilson||T999| [Milam et ar||2005| . 



These isotopic abundances are typical of what is measured in nearby starbursts (eg. Henkel et 
al. 1994), and are within a factor of ~2 of the entire range observed in the Galaxy. A ^^C0(1- 
0)/^^CO(1-0) line ratio of 3.0 implies ti^co{i-o) — 1) these isotopic abundance ratios. 

Uncertainties in derived opacities depend sensitively on the true [^^00/0^*^0] abundance ratio 
which may differ from the value adopted here. Both [CO/^^CO] and the [CO/C^^O] decrease with 
stellar processing (assuming the CO isotopologues abundances are proportional to their respective 
isotopic abundances). Based on Galactic disk studies, [^^CO/C^^O] is expected to decrease from 
~7.5 in the local ISM to ~4 in the inner kpc of the disk, arguing for a decrease in the combined 



ratio with nuclear synthesis (eg. Wilson 1999; Milam et al. 20051. Attempts to determine isotopic 
abundances directly in starbursts obtain ^^CO/C^^O ~ 5, further supporting a lower ratio in 



high metallicity regimes (eg. Henkel et al. 1994 1 . On the other hand. Galactic Center (Sgr B2) 
determinations arrive at anomalously high values of [^^CO/C^^O] ~10 (eg. Langer &: Penzias|l990 1. 
Using the Galaxy as a guide [^'^CO/C^^O] should fall somewhere between 4 and 10. We favor values 
on the low end of this range for the high metallicity nucleus of Maffei 2 for three reasons: (1) The 
Galactic disk gradient and the starburst values imply low ratios in heavily processed locations. (2) 
The LVG models give consistent ^^CO and C^^O parameter space solutions for value [^^CO/C^^O] 
~ 4, but not for 10 (Q. (3) There is marginally significant evidence for lower CO(1-0)/^^CO(1-0) 
and ^^CO(1-0)/^®CO(1-0) line ratios along the central ring, even towards the lower column density 
portions (Figure [g]). This may represent direct evidence for enrichment of C^^O relative to ^^CO 
(and CO) in the immediate vicinity of the nuclear starburst (similar effects are seen in IC 342; 
Meier Turner[[2001[ ). 



Modulo small differences in resolved-out flux or linewidth, we conclude Ti-^(^o(i-o) ~ 1 over the 
molecular peaks, but note that systematic abundance uncertainties allow anything between r <C 1 
up to r ~ 4. The fact that the CO(3-2)/^3(2!o(i-0) line ratios of [Hurt et ari ( [l993b[ ) imply much 
lower opacities must then be a result of CO(3-2) emission (and likely CO in general) preferentially 
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originating from higher excitation gas than do ^^CO(l-O) and ^^CO(l-O). 

With Ti3co(i-o) independently constrained from the -'^^CO/C-'^^O ratio, comparisons between 
Tex derived from eq. ([l]) (corrected for resolved out flux) and the observed ^•^CO(l-O) peak bright- 
ness temperature, ^^Tmb, gives constraints on the filling factor, ^^/, via ~ [^^Tmf,/^^ Jio(7ea)](l — 
g-n3co(i-o)^-i_ Towards the molecular peaks ^'^/ ~ 0.33 is estimated. Given the potentially large 
uncertainty in the estimate of the average ^'^CO opacity, ^^f should be considered only indicative. 
If ''"i3CO(i-o) ^ 1 then could be as low as ~0.15. The relatively bright ^^CO(l-O) emission and 
the fact that ^^f < 1 requires that ti^(^o{i-o)^ 0.33. 

In summary, the J= 2-1 to 1-0 line ratios of the lower opacity CO isotopomers imply LTE 
excitation temperatures of Tex ~3-10 K. Brightnesses of the higher opacity J > 1 transitions of 
CO, suggest that they preferentially sample more limited volumes of warmer gas. The opacity of 
the ^^CO(l-O) transition appears to approaches unity over much of the nuclear peaks. 



5.3. HCN 



We also compare the distribution of dense gas traced by HCN to that of the total gas traced 
by CO. Figure [5|:, shows the CO(1-0)/HCN(1-0) integrated intensity line ratio. The ratio varies 
from ~8 at CMC F to >20 at the ends of the molecular bar. Figure ^ shows the general radial 
trend commonly seen in galaxies, namely an increase in the CO/HCN intensity ratio as one moves 



away from the star forming sites at the center (e.g., [Heifer Sz Blitz 1993, 1997[ Sorai et al. 2002 1 . 
Sites of high CO(1-0)/^'^CO(1-0) ratios (particularly in the off-arm regions) also have the highest 
CO(1-0)/HCN(1~0) ratios (Tablej?]). This provides evidence that the ^^00(1-0) is at least partially 
sensitive to the density of the gas. 



6. Nuclear Gas Kinematics: The Parallelogram and a Bar Model 



Maffei 2 is a strongly barred galaxy with a disturbed morphology, probably due to interaction 



with a nearby companion (Hurt et al. 1993a Hurt, Turner &; Ho 1996; Mason &; Wilson 2004 1 . 



Position- Velocity (P-V) diagrams based on lower resolution CO(l-O) data, show that the molecular 
gas in the central regions reaches high (^J 75 km s~^) radial velocities over very small projected 
radii ( [Ishiguro et"aL]|1989[ [Hurt Turner ||1 99 1| ). |lshi guro et al. ( |1989[ ) has interpreted this feature 
as an expanding ring generated by an explosive event some ~ 5 x 10^ yrs ago, superimposed on a 
Keplerian component. Since Maffei 2 is strongly barred, it is reasonable to consider whether these 



motions are a result of non-circular motions due to a barred potential. Ishiguro et al. (19891 argued 



against this based on the fact the position angle of the "molecular bar" and the major axis of the 
galaxy are close so that any non-circular motions from a bar would be in the plane of the sky and 
therefore could not explain the motions observed. However, if there are ILRs (so X2 orbits exist) 
and / or a nuclear bar rotated with respect to the large scale bar exists then this is not the case. 
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The gas kinematics can been seen in the channel maps of ^^CO(l-O) (Figures g and g and 
HCN(l-O) 1^00(1-0) traces the velocity field of the total column density, while HCN the 

velocity field of the dense gas. CO emission extends from vlsr ~ —160 to +100 km s~^ with 
blueshifted emission in the north. HCN is confined to velocity ranges vlsr ~ —145 to +60 km 
s~^. If we assume trailing spiral arms, the northern arm is the near arm, consistent with the larger 
internal extinction there (Figure [2]: 



Position-velocity (P-V) diagrams for Maffei 2, shown in Figure 10 were made from the cubes by 
averaging the central 5 along the major axis of the galaxy {p.a. = 206°; Table[l|. We constructed 
major-axis P-V diagrams for CO(l-O), 13CO(1-0) and HCN(l-O). The CO(l-O) and ^^CO{1-0) 
P-V diagrams reveal an essentially complete "parallelogram" associated with the central ring. The 
parallelogram has a width of ~10", a spatial extent of 20"(160 pc in radius), and a total velocity 
extent of ~250 km s~^ (uncorrected for inclination). The ^^CO-emitting gas along the central ring 
appears to be nearly uniform. CO(l-O), on the other hand, is asymmetric with much brighter 
emission at the CMC D + E starburst sites along the eastern side of the ring. The cause of 
the asymmetry is unclear, but since CO is so optically thick, its emissivity is more susceptible to 
locally elevated kinetic temperatures or other non-LTE surface effects. The uniformity of the ^^CO 
is likely to be a better indicator of gas surface density in this situation. Beyond the parallelogram, 
the velocity field is dominated by two peaks corresponding to the ends of the molecular arms. 

The velocity field of HCN(l-O) is significantly different from that of CO. The parallelogram is 
not apparent in the HCN P-V diagram, but instead dominated by two peaks corresponding to the 
intersection of the molecular arm emission and the parallelogram in the P-V diagrams. Even along 
the central ring, HCN has a much lower covering fraction in velocity space than CO. 



6.1. Bar Model for Maffei 2 



The new CO P-V diagrams (Figure 10) have sufficient spatial resolution to reveal that the 



"oval"-shaped pattern is actually a "parallelogram" like that observed towards the Galactic Center 



(Bally et al. 1988 Binney et al. 1991). The similarity of molecular gas kinematics in Maffei 2 to 



the Galactic Center, which is explained by gas response to a barred potential (eg. Binney et al. 



1991 Huettemeister et al. 1998 Rodriguez- Fernandez et al. 2006), leads us to construct a model 



of barred gas response in Maffei 2. 

We have modeled the gas distribution and kinematics in response to a stellar bar using an 
analytic weak-bar model. Such models are based on treating the dissipational nature of gas with 



the addition of a damping term proportional to the deviation from circular velocity (Wada 1994 



Lindblad & Lindblad 1994; Sakamoto et al. 1999). Despite their simplicity, these model matches 



the structures seen in full hydrodynamical simulations with surprisingly fidelity (eg., Lindblad Sz 



Lindblad|1994 ). The simplicity of an analytic bar model permits us to quickly explore a wide range 
of bar parameters. 
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Our model is based on those of Wada (1994) and Sakamoto et al. (1999), with the following 



modifications. (1) The gas dissipation term is extended to include azimuthal damping. This is 
done by adding a term, 2fi(j)i, to eq. (2) of Wada (1994) analogous to their eq. (1). (2) The model 



is extended to include both a small scale (hereafter "nuclear") bar and a large scale bar, according 
to the prescription of Maciejewski &; Sparke ( |1997 2000). (3) The axisymmetric potential has been 
changed so that it generates a double Brandt rotation curve coupled in size to the two bars. The 
second component of the potential is required both to match the flattening of the rotation curve 
near the nucleus ( Hurt, Turner Ho|l996 | and to generate the necessary existence of ILRs at ~5". 
The size and shape of the large scale bar is set to match the NIR, CO and HI images of the galaxy 



(Figure 11 1; Hurt et al. 1993a: Mason k Wilson 2004 Hurt, Turner k Ho 1996). It is assumed 



that the nuclear bar is co-planar with the large scale disk. 



Similarities between the model and both the true gas distribution (Figure 1 1 ) and kinematics 
(Figure 12 &; 13 ) are excellent. Table[9]lists the fitted model parameters. In the weak bar scheme the 
molecular gas is predicted to follow the sites of orbit crowding (higher density of dots in the figure) . 
The model velocities match the observed pattern quite well. In general, the models are robust to 
small changes in parameters as long as two main requirements are met, (1) the combination of 
potential/rotation curve and bar parameters are such that there are two nuclear inner Lindblad 
resonances, oILR and ilLR, and (2) the size scale of each bar is close to the values chosen. The 
first is important because it is this condition that is required for perpendicular X2 orbits. In 
our barred model, the perpendicular X2 orbits are vital to explain the large 1-o-s velocities and 
"parallelogram" feature seen close to the center. The second point is important because it sets the 
scale of the features seen in the gas (ie. the xi orbits run between the oILR and corotation). 

Our model shows gas associated with the well known xi and inner perpendicular X2 orbits of 
barred potentials (eg. Athanassoula|[l992 ) . A good fit for the nuclear morphology is achieved when 
the 4:1 ultraharmonic resonance of the nuclear bar is set equal to the oILR of the large scale bar. 
The pattern speed of the nuclear bar in Maffei 2 is thus much higher than the pattern speed of the 
larger bar, implying that the nuclear bar is decoupled from the larger scale bar. The morphology 
is best matched with the nuclear bar rotated ~ 10° clockwise (as viewed from the perspective of 
Figure [TTfe) relative to the large scale bar. 



The "parallelogram" of the ^^CO P-V diagram reveals additional information about the gas 
associated with the closed nuclear bar orbits. The ^^CO parallelogram is nearly complete. This 
suggests that, while not obvious in the integrated intensity map due to the high inclination (though 
evidence is seen for it ~ 5" west of GMC E), the entire oval X2 orbit region appears to contain 
molecular gas. The emission along the central ring in the P-V diagram is nearly unresolved spatially, 
suggesting that only a very small range of X2 orbits are populated. Contrasts of >6 in column 
density are seen between the molecular gas associated with nuclear X2 orbits and the very center 
of the galaxy, implying that the vast majority of inflowing molecular gas is trapped at the central 
ring and does not reach the core of the galaxy. The parallelogram suggests that a majority (~60 
%) of the ^^CO emission originates from gas residing on X2 orbits with the rest residing on the xi 
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orbits. 



The observed "parallelogram" is slightly wider along the position axis than the model predicts. 
This suggests that the inner X2 orbits are slightly more circular than displayed in Figure [TT]i. This 
is likely due to the limits of the epicyclic approximation at the very center of the potential. Nearly 
circular X2 orbits are a common feature of the more complete hydrodynamical simulations (e.g., 



Finer, Stone k Teuben 19951. 



7. Molecular Clouds in the Nuclear Environment: Cloud Properties as a Function 

of Location 

Armed with a reasonable kinematical model of the center of Maffei 2, and having identified 
the sites of current star formation from mm continuum maps, we can investigate the effects of 
environment on the properties of molecular clouds. Densities, n/^j, and kinetic temperatures, 
Tfc, of the individual nuclear GMCs were determined by running Large Velocity Gradient (LVG) 



radiative transfer models with the observed intensities and line ratios as inputs (eg. Goldreich &; 
Kwan||1974[ [Scoville Solomon|[T974} |De Jong et al]|1975[ ) . 



We adopt single component LVG models for the clouds. Three independent parameters, 
nH2i Tfc and Xco/dv/dr are varied over the ranges n/fj = 10^-10^ cm~^, = 1.5-150 K, and 
Xco / dv / dr = 10"'^ — 10~'''^ (collision coefficients are from De Jong et al.|[l"975 ). For each location. 



seven (eight if including HCN) distinct measurements, the two isotopic line ratios, the two A J line 
ratios, the peak T^fe of the uniformly weighted CO(l-O) map and the ratio of cloud linewidth to 
cloud size are used to constrain the model parameters. The itlcj ranges do not include system- 
atic uncertainties associated with changes in Xco/dv/dr or more general uncertainties related to 
the validity of the LVG approximation itself. Additional model solutions with Xco/dv/dr varied 
by a factor of ±0.3 dex were determined (not shown). Increasing the velocity gradient (or corre- 
spondingly decreasing the abundance) results in an increase in derived densities by ~0.3 dex and a 
decrease in derived T/^ of ~ 10 K. This variation is indicative of the sensitivity of the physical con- 
ditions to changes in the abundance or velocity gradient at the level permitted by their systematic 
uncertainties 



In Figure [T4| LVG model solutions for six locations across the central bar are displayed. Typ- 
ical values of the velocity gradient for the GMCs are ~ 1-2 km s"^ pc"^ (Table |4|. We force 
abundance per velocity gradients of, Xco/dv/dr = 10~^'^^ (10"^ ''^) for ^^CO (C^^O), correspond- 
ing to [CO/H2] ~ 8 X 10-^ [C0]/p3C0] ~ 60, [C0]/[Ci80] ~ 250 (fg, and dv/dr ~ 1.5 km s'^ 
pc^^ (Table [4]). The itlo" range for each ratio and the measured value for the CO(l-O) antenna 
temperature constrain parameter space. 

Average densities and kinetic temperatures implied by the LVG models are ~ ]^q2. 6-3.0 
cm~'^ and T^ ~ 15-35 K. All the mapped CO lines imply a consistent set of physical conditions. 
Densities derived from CO tend to be nearly constant across the nuclear bar. Clouds associated with 
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the starbursts (D and E) are slightly warmer than the others {Tk ~ 30-40 K). By contrast, slightly 
cooler and denser values are derived for the quiescent gas clouds on the nuclear X2 orbit ('Western 



Ring'; Figure 14) compared to the starbursting side of the central ring (GMC E). The solutions 
reproduce the observed brightness temperature of the uniformly weighted CO(l-O), indicating a 
filling factor of unity. A unity filling factor for CO(l-O) is consistent with that estimated from ^^CO 
excitation given its larger beam size. Only GMCs D and F have predicted brightness temperatures 
slightly lower than observed (by ~50 %), which could be explained if the CO-emitting surfaces are 
somewhat warmer than the bulk of the gas in these clouds. The ^^CO and C^^^O solutions are in 
excellent agreement, suggesting that the adopted relative abundances are reasonable. 



HCN LVG models were also run, for levels up to J=12 (collision coefficients from Green k. 



Thaddeus||l974[ ). Overlaid on Figure [H] are the observed range (ilcr) for the ^3CO(1-0)/HCN(1-0) 



line ratio, derived from the HCN(l-O) models. An abundance per velocity gradient, XijcTv/dv/dr 
= 2 X 10~^ was assumed, consistent with Galactic HCN abundances and dv/dr = 1.5 km s~^ (eg., 
Irvine, Goldsmith &: Hjalmarson|1987 Paglione et al.|1998 1 . A small correction for resolved out flux 



has been made assuming emission is uniformly extended on scales i^30 . Densities derived from the 
HCN(l-O) line are about an order of magnitude higher than fit from the CO isotopologues: HCN 
is brighter than expected based on the CO-derived physical conditions. This implies that (1) these 
molecular clouds have a significant component of denser clumps from which the HCN originates, or 
(2) the HCN abundance is much larger than has been adopted. The absolute abundance of HCN 
is not known well enough to eliminate the second possibility, but the magnitude of the increase 
required (i^lOx) makes it unlikely. The HCN data suggests that while the single-component LVG 
approximation yields internally consistent solutions for the optically thinner isotopologues, which 
sample the bulk of the molecular column density, it breaks down when including the very high 
density gas. The derived nj^j ^-nd T^ should then be treated as volume averages for the gas clouds 
traced in the isotopologues, but not necessarily the whole of the ISM. The ^^CO/HCN would then 



reflect the relative fraction of very dense gas at each position (e.g., Kohno et al. 1999[ Meier &; 
Turner||2004l ). 



To summarize, molecular clouds in the central 300 pc of Maffei 2 averaged over ~60 pc scales 
tend to be only modestly denser than GMCs in the disk of the Galaxy. How much denser depends 
on the exact velocity gradient /abundance of CO and HCN present. If the molecular gas has a 
velocity gradient ~ 1-2 km s~^ pc~^, consistent temperatures and densities are obtained from the 
CO emitting gas, at values of < n/^j > — 10^'''^ cm~^ and T^ ~ 15 - 35 K for all GMCs. At these 
densities, the CO isotopologues are subthermally excited. The HCN emission implies subclumping. 



7.1. CO as a Tracer of Molecular Gas Mass: The Conversion Factor in Maffei 2 

There are indications that CO(l-O) is overluminous per unit mass of H2 gas in the nuclear 
regions of spiral galaxies relative to Galactic GMCs, and thus use of the Galactic conversion Xco 
can lead to overestimates of molecular gas masses in these systems, (e.g. Dahmen et al.|1998 Meier 
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Turner|200T| [Weiss et al.|2001t [Meier fc Turner |2004[ ) . In this section, we compare several different 
methods of estimating molecular gas column densities to assess the validity of the conversion factor 
in the nucleus of Maffei 2. 



7.1.1. Molecular Gas Column Densities from the Optically thin CO Isotopologues, Dust 

Continuum, and the Virial Theorem 

Optically thin lines of CO isotopologues allow estimates of the molecular gas column density 
directly by summing up the emission from each molecule. These estimates depend only on the 
knowledge of relative CO abundance and excitation. For LTE (eg. Scoville et al.|[1986 |, 



, , 17 o iTex + 0.921 \ H2] 

N{H2)^co = 5.75 X 10^^ cm~^ ^ ^ |^ e-e. 



^ , ^^co {K km s-'). (2) 

where p^^jj is the abundance of the isotopologue, V^, V and *To are the frequency (in GHz), 
the opacity and the characteristic temperature {hvo/k) of the particular transition (*To = 5.29 for 
^3CO(1-0) and 5.27 for 0^*^0(1-0)). We calculate Tex separately from the ^3cO(2-1)/13qq(^.o) 
line ratio assuming that ^^tco — 1 (^) and the C^^O(2-1)/C^^O(1-0) line ratio. These values 
and the H2 column densities derived from them are presented in Table [s] for each peak (Tex = 5 K 
is assumed for positions beyond the (2-1) primary beam). Under LTE, T^^x values turn out to be 
almost independent of opacity over the range of - 5 (changing by < 20 %) when J = 2-l/J=l-0 
line ratios are around the observed value of ~ 0.7. Therefore systematic uncertainties in stem 
primarily from uncertainties in the assumed isotopic abundances. Abundances are expected to be 
within a factor of 1^2 of the adopted values (^. 

Values of N^a range from < 0.61 (1.3)-10 (6.5) x 10^^ cm'^ (Tableg, based on ^^CO (C^^O) 
fluxes, with corresponding mass surface densities of S ~ <130 (280)-2200 (1400) Mq pc~^. 
Emission from the fainter C^^O line does not extend to the lower column densities, otherwise the 
predictions of column densities from the two species agree to within a factor of 2. 

Dust continuum emission has also been detected at A =1.4 mm towards the several of the 
GMCs, which gives another estimate of the molecular gas mass. After accounting for the free-free 
contribution (Table [s]), dust fluxes, SiAmm^ range from 7-19 ±3 mJy for each cloud. Assuming a 
gas to dust ratio of 100 by mass, the gas mass is related to the 1.4 mm dust continuum flux by 



(eg., [Hildebrand[[l983[ ): 

AWl.^.-) = 300Me(%|f) (^)^ 

(3) 



cw? g 1 



- 10.50 
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where is the dust absorption coefficient at this frequency, Si,4^mm is the 1.4 mm dust flux, D 
is the distance and is the dust temperature. The dust opacity, k^,, at 1.4 mm is taken to be 
3.1 X 10~^ cm^ g~^, uncertain by an estimated factor of four (Pollack et al. 1994). We adopt 
Trf = 40 K based on the FIR colors ( Rickard &: Harveyl|1983 1. The dust temperature applicable to 



the 1.4 mm observations could be lower than this if a cool dust component undetectable shortward 
of 160//m exists. The existence of a cooler dust component would cause us to underestimate the 
implied molecular gas mass. This is likely more important away from the nuclear region. Dust 
masses for the cloud peaks are listed in Table Isj 



Virial masses can be derived from linewidths, following the treatment of Meier Sz Turner (2001 1 



since the individual GMCs in the center of Maffei 2 are resolved. Virial masses are given in Table 
[8j They have an intrinsic uncertainty of about a factor of two due to internal cloud structure. In 
addition, if systematic motions such as cloud streaming motions or non-circular bar motion, are 
present within a single beam (almost certainly the case; the linewidths due to the internal 
gravity will be overestimated. In short, the virial masses will be upper limits to the true cloud 
masses. 

Finally, as a crude consistency check, a column density, ^Nh2-, is calculated from the LVG 
model derived densities, by averaging the number density over an assumed depth of \fW^h- These 
values are also recorded in Table [Sj These values represent upper limits to N(Il2) if nn^ is confined 
to a fraction of this volume. 



7.1.2. The CO Conversion Factor in Maffei 2 

We can compare the three different column density estimates — optically thick CO {^'^"Nh^), 
optically thin ^^CO (^^Nh^) and C^®0 (^^iV/fa), and dust {'^Nh^) — to estimate a CO conversion 
factor, Xcoi for the nucleus of Maffei 2. Column densities based on the CO isotopologues (Table |8j) 
are lower than the those derived from CO(l-O) intensities using the Galactic value of Xco — 2 x 10^*^ 
cm"2 ("k: km s"1)-^ ( [Strong et al.|l988[ [Hunter et al.|l997{|Dame et al.|200l| ). Values from the thin 



C^^O lines are ~2-4 times lower than the Xco estimates. If we require that II2 column densities 
derived from opacity-corrected ^^CO(l-O) and C^^O(l-O) agree (thereby constraining Ti^co{i-q)) 
then Galactic values of the conversion factor can be reached only for [CO/C^^O] > 600. Given the 
high metallicity environment of the nucleus of Maffei 2, this seems unlikely. Away from the nucleus 
there is some evidence for another factor of two further decrease in the conversion factor; however, 
statistical uncertainties are at least this large due to weak emission and Tg^^ not being determined 
towards these locations. 

Uncertainties estimated for the gas column derived from dust emission are higher than for the 
isotopologues, but they too tend to support lower gas columns than predicted by the Galactic Xco- 
Dust-based gas masses are also lower than the Xco values by factors of ~2-4 for the adopted dust 
parameters towards the detected GMCs (Table [8|). Gas column densities estimated from the dust 
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are in good agreement with the opacity corrected ^^CO estimates except for GMC E, but trend 
a factor of ~30 % higher than those from the C^^O isotopologues. This is an indication that the 
uncertainty in these column densities are at least this large. While these methods are different, we 
do not claim they are completely independent, because there may be hidden correlations between, 
say, CO relative abundance and dust to gas ratio. But the dependences on metallicity and other 
factors such as temperature are not necessarily the same for these mass tracers. That the gas 
column densities estimated from the dust and C^^^O are both low provides additional confidence 
for the assertion that the gas column densities are overestimated by the Galactic value of Xco- 

Column densities obtained by averaging the virially derived masses are higher than the other 
methods, which is not surprising. Linewidths in the central region (particularly GMCs C, D, E 



and F) include two distinct components moving on completely different orbits (Figure 10), and 
so systematic motions within the ~60 pc (line-of-sight) beam due to the bar orbits ( ^6.1[ ) cause 
an observed linewidth larger than random gravitational motions within the clouds would imply. 
Because of the presence of this motion, we expect that virial methods using observed linewidths to 
be severe overestimates of the cloud masses. 

In summary, we conclude that the conversion between ^"^00(1-0) and H2 column density ap- 
plicable to the central region of Maffei 2 is ^'^'^Xco —0.5 -1.0 xlO^° cm~^ [K km s~^)~^ , ^2-4 
times lower than the Galactic value, with uncertainties of ^100 %. 



8. A Nuclear Bar-Driven Starburst in Maffei 2 

8.1. Star Formation Rates and Efficiencies 

The rotation curve from the double bar model can be used to estimate dynamical masses 
directly without having to try to remove the non-circular velocity component (Figure [Tsja). The 
dynamical mass over the central ring is y^dyn{R = 7") = 2.1 X 10^ M0. The molecular mass 
estimated from C^^O over the same region is 6.9x10^ Mq. Dynamical masses for the central 20 
radius are Mrfj^„(20") = 7.3 x 10® -^0, while the molecular mass over this region is 2.1x10'^ -^0- 
Molecular mass fractions are thus ~3% percent over much of the central molecular bar. Molecular 
mass fractions scale as the distance, so the uncertainty in the distance to this galaxy (§2) can 
change these values by up to a factor of 2. Resonant structures, such as the molecular bar observed 
in the nucleus, are probably driven by the stellar potential rather than the gas. 



Lyman continuum ionization rates, ^lyc (for Te = 10 K; e.g. Mezger &; Henderson 19671 



and star formation rates based on the 89 GHz continuum are given in Table [6| To produce the 



total observed free- free emission across the central 30 , the excitation of 2600 effective 07 (Vacca, 



Garmany &; ShuU 1996) stars is required. A significant fraction of this ionizing flux (~ 1.4x 10^^ 

or ~ 1400 "effective" 07 stars) arises near the two central molecular peaks (GMCs Dl-I-E and F). 

Towards radio continuum sources, I, II and IV, the local star formation rates are 0.05, 0.05 and 
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0.04 Mq yr ^, respectively, based on the conversion between N^yc and SFR of Kennicutt (19981 



These values match the star formation rate predicted from the HCN(l-O) luminosity using the 



relationship that Gao k, Solomon (2004 1 have derived from large scale HCN measurements. The 



relationship between HCN(l-O) luminosity and star formation on 60 pc scales in Maffei 2 is the 
same as that observed on kpc scales in luminous infrared galaxies. 

The ionization rate across the nuclear bar corresponds to a massive star formation rate of ~ 
0.26 Mq yr-i with ~0.14 M© yr ^ originating from the nuclear ring. At this rate the molecular gas 
in the ring could sustain the current star formation rate for ~ 5 x 10'' yrs if no gas replenishment 
from the arms occur. If a ZAMS Salpeter IMF with an upper (lower) mass cutoff of 100 Mq (0.1 
Mq) is adopted then a total stellar mass over the central ring of M^, = 5.4 x 10^ Mq is generated in 
the current burst. These values correspond to star formation efficiencies, SFE = M^I{M^ -^Mh^^ 
of ~10% over the central ring, peaking at the nuclear x\ - X2 orbit intersections. The SFE drops 
to ^ 4% along the molecular arms, similar to Galactic disk values. 



8.2. Gas Inflow, Stability and Triggered Star Formation in Maff'ei 2 



What drives the star formation in the nucleus? Is it the large molecular gas surface density 
or is there evidence for a trigger that is unique to the nuclear region? What is responsible for the 
large concentration of nuclear gas in Maffei 2? With these observations we can address the link 
between star formation and molecular gas on GMC sizescales in the nuclear region of Maffei 2. 

In the context of our dynamical model, the presence of the large gas mass is probably due 



to slow inflow along the nuclear bar (eg. Roberts et al. 1979 Athanassoula 1992; Turner & Hurt 



1992 Regan, Vogel, & Teuben 19971 iSheth et al. 20051. Is there sufficient gas inflow to produce 



the observed star formation? It is assumed that the inflowing gas will form stars, and that the 
star formation process is initiated at the location of the xi-X2 orbit intersection in the nuclear ring. 
Then the radial gas mass flux at this galactocentric radius determines the star formation rate. The 
gas mass flux is related to the average inflow velocity, Vinf, the average arm mass surface density, 
Tiarm, and the arm width, w. From the ^'^CO data Harm — 210 Mq pc~^ and w ~ 5" (80 pc). 
Inflow velocities are determined from the bar model. Typical values are -20 to -40 km s~^ along the 
bar arms. Averaged over the arm area only, Vr — -21 km s~^. Adopting these values (an upper 
limit), a mass inflow rate, Mj„j ~ 0.7 Mq yr^"*^ is obtained. Since Mj„j is a factor of ~5 larger 
than the nuclear ring star formation rate estimated from the millimeter continuum, the inflow rate 
is sufficient to fuel the nuclear starburst even with modest efficiency. 

Does the molecular gas form stars due to gravitational instabilities or is it directly triggered? 
The gravitational stability of a thin, rotating disk can be estimated from the Toomre Q parameter 

< 



( Safronov 



1960 



Toomre 1964 1 . A gas disk is unstable to gravitational collapse if Q = kct/ ttGT, 



gas 



1, where k is the epicyclic frequency, a is the gas velocity dispersion and S^as is the gas surface 
density. Figure [Tsb displays the azimuthally averaged values of k from the bar model together 
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with the observed ^^CO velocity dispersion, mass surface density, and corresponding Q values. Q 
is 8-10 across the central ring region containing the starburst and remains >1 over the central 30" 
radius. Clearly the data are not consistent with star formation occurring in Q < 1 gravitationally 
unstable gas. This is not surprising for gas in the very center of galaxies given the (1) strong 
noncircular motions present, (2) the failure of the thin differentially rotating disk approximation 



and (3) potentially strong turbulence and magnetic fields (e.g., Elmegreen 1999 Combes 2001 



Wong & Blitz 2002). 



Another estimator for gravitational stability that may be more suitable to nuclear gas can be 



obtained from Elmegreen (1994|. Elmegreen (1994) estimates the critical density above which gas 
in a ring associated with an ILR can collapse to form stars as pcnt = 0.6 k^/G, or ricrit = 2.08 x 
10"^ K^(km s~^ kpc~^). The epicyclic frequency at the radius of the ring is k ~ 2000 km s~^ 
kpc~^ which implies ring densities must be n}{2 > 8 x 10^ cm~^ to be unstable to collapse. From 
the LVG analysis we find that the average density of the CO-emitting gas along the central ring is 
an order of magnitude lower than this value. 

A lower limit to the stability of the molecular clouds can be set by assuming the clouds 
remain gravitationally bound against tidal forces. A cloud of mass, m, will remain bound if 



m 



M{R){r/Rf, where M{R) is the total mass enclosed within a galactocentric radius, i?, and r is 



the size of the molecular cloud (eg.,|Stark et al. 1991 ). Clouds with densities oin^^"-^ ^ 3.6Mmq{R)R, 



"pc 



remain bound. For R = 80 pc and Mjv/q(80 pc) = 1.1 x 10^ Mq, values applicable to Maffei 2's 
nuclear ring, n''^'' ~ 630 cm~'^. This value is close to the densities inferred from our LVG analysis. 
The average molecular gas densities in the central ring are too low to be gravitationally unstable, 
and are likely only marginally tidally bound. Indeed, along much of the central molecular ring not 
associated with the sites where the arms terminate, little star formation is observed. 

It appears, then, that gravitational instability is not the answer. Instead we consider the 
possibility that the star formation is triggered by events external to the clouds. Star formation in 
the nucleus of Maffei 2 is concentrated at the location of the xi-X2 orbit intersections indicated by 
our modeling. At these xi-X2 orbit intersection regions star formation appears to be triggered by 
the collision of gas flowing inward along the arms of the bar with the existing, more diffuse gas of 
the central ring. 

We propose that the evolution of the nuclear starburst has proceeded as follows. A recent 
interaction between a small companion and Maffei 2 has driven a large quantity of gas into the 



nucleus, building up a compact central bulge seen in the NIR (Hurt et al. 1993a; Hurt, Turner Sz 



Ho 1996). Assuming the potential generated by this compact bulge is slightly oval (a few percent 



is all that is necessary), it has forced the nuclear molecular gas into the bar distribution currently 
observed. Inflow along the nuclear xi orbits piles up gas at the nuclear xi — X2 orbit intersections. 
The interaction results in a fraction of the molecular gas going to the formation of dense cloud cores 
which collapse and trigger the star formation events at CMC D and just downstream of CMC F. 
Gas not incorporated into the dense component at these locations is then tidally sheared into the 
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moderate density, nearly uniformly distributed ring, which in turn becomes the target for future 
collisions with infalling gas. As long as there is gas flowing inward the burst of star formation at 
the arm-ring intersection can be sustained. 

This scenario provides a good framework for all of the molecular gas and millimeter continuum 
observed toward Maffei 2, with one exception, the star formation near GMC E. This star- forming 
complex is on the X2 orbit but not at either of the X1-X2 intersection regions. Nor is it a strong 
HON source (though some HON emission is seen). Why is star formation occurring here? Two 
possibilities come to mind: (1) The star formation is triggered by the molecular gas associated 
with GMC E interacting with the molecular gas towards GMC D after having traversed one half 
of the X2 orbit. (2) The star formation here reflects a slightly earlier epoch event associated with 
its passing through the southern xi - X2 interaction region. It is now being seen with a time lag 
equal to the traversed portion of the ring divided by the orbital velocity. From the nuclear ring 
parameters the time lag would be ~ 1 Myr. That the spectral index of the millimeter continuum is 
somewhat steeper towards GMC E, possibly suggesting a slightly older starburst with more evolved 
and less dense H II regions, may favor the latter scenario. 



8.3. Comparisons of Maffei 2 with Other Nearby Nuclei 

The nuclear morphology of Maffei 2 is similar to that observed in the bright nuclei of the barred 



Regan & Vogel 


1995 


Sakamoto 


Binney et al. 


1991 


] 


ilodriguez- 



Fernandez et al. 2006 ) . All have nuclear bar morphologies reminiscent of their large scale analogs 
(e.g., |Athanassouia|[T992l ). In IC 342 and M 83, it remains somewhat ambiguous whether they are 
nuclear bars or just the inner portions of the large scale bar, due to a combination of having massive 
clusters that potentially influence the dynamics (e.g., Schinnerer et al.|2003 Crosthwaite et al.|2004 



Sakamoto et al.|2004 Schinnerer et al.|2007 | and low inclination, which hampers kinematic studies. 
Maffei 2's kinematics leave little doubt that it is a true double bar. In fact the CO velocity field 
in the nucleus of Maffei 2 is perhaps the best current example of nuclear non-circular, bar motions 
outside our own Galactic Center. Therefore Maffei 2 can be added to NGC 6946 as confirmed 
double barred galaxies, but with a physical scale about three times larger. It is interesting that like 
NGC 6946 (and NGC 2974; [Ishizuki et aI][T990l [Emsellem et~al][2003l [Schinnerer eTal][2006| , our 
CO(l-O) observations imply the existence of straight shocks in nuclear bars. The inner ring and 
offset straight shocks do not appear to be common features of hydrodynamical models of secondary 



bars (eg., Shlosman & Heller 2002 Maciejewski et al. 2002 1 . Moreover, this conclusion seems to 
hold true for galaxies with both strong large scale bars (Maffei 2) and weak large scale bars (eg. 
NGC 6946). 

These nuclear bars also influence physical conditions of the gas. Despite the presence of 
luminous star-forming complexes in these nuclei, emission from the lines of the CO isotopomers is 
dominated by subthermally excited emission from low excitation (Tg^ ~ 5 — 15 K) gas, and that 
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this emission represents the properties of the bulk of the molecular gas. While the ISM in Maffei 
2 appears slightly warmer than IC 342 (Meier &; Turner 20011 its are very close to the average 



properties of NGC 6946 and the outer gas lanes of the Galactic Center ( Huettemeister et al.||1993 



Paglione et al.|[l998 Meier &: Turner||2004 Nagai et al.||2007 |. However, densities in the nucleus of 



Maffei 2 are consistently about 0.5 dex lower than NGC 6946 and the Galactic Center. We suggest 
this comes from Maffei 2 having a stronger nuclear bar than the other two nuclei, resulting in more 
dramatic disruption and redistribution of its nuclear ISM. 



9. Summary 

New aperture synthesis maps are presented for emission in the J=2-l and 1-0 transitions of 
^■^CO and C^^O, as well as the J=l-0 lines of HCN and CO in the central arcminute (~ 1 kpc) 
of Maffei 2. The II2 column density as traced by optically thin CO isotopologues is similar in 
morphology to what is implied from ^■^CO, except that the emission from the isotopologues is more 
closely confined to the two extended molecular arm ridges and more uniformly distributed across 
the central ring. The dense gas traced by IICN(l-O) is more confined to the center of the galaxy 
than the CO emitting gas. 

The central molecular bar contains five main peaks that resolve into at least 17 distinct GMCs, 
with radii of ~40-110 pc and linewidths ^ 40 km s^^. In the two innermost molecular cloud 
complexes, at galactocentric radii of ~ 5" (80 pc from the dynamical center), the GMCs are 
distinctly nonspherical, elongated along the nuclear bar, with linewidths as large as 100 km s~^. 
These GMCs are probably being tidally stretched due to the nuclear potential. 

The H2 column density for the central GMCs is ^ 4.4 - 10 x 10^^ cm'^ (S ~ 950 - 2200 
M0 pc~^), corresponding to mean optical extinctions of Ay ~40-100. The molecular mass within 
the central 20" galactocentric radius (~300 pc) is 2.1x10^ Mq, while the dynamical mass in the 
same region is M(^y„(20") = 7.3 x 10^ Mq. The molecular mass is only a few percent of the 
dynamical mass. Excitation temperatures, assuming ^^r(^^r) ^ 1 (^1), 3,rG ^ ex 

3 - 6 K 

over much of the central 500 pc for both ^'^CO and C^^O. These T^x values are low compared with 
the brightness temperature observed in CO (s;30 K) indicating subthermal excitation, and that the 
average densities of the GMCs are probably only moderate. Single component LVG analysis of the 
GMCs in CO, ^^CO, and C^^O yield best-fit solutions of n//^ ~ lO^-^s cm'^ and Tkin ^ 20-30 
K. Average densities estimated from the total C^^O column densities are consistent with these 
values. 

The ^^CO and C^^O lines are weaker than expected from CO(l-O), which appears to be over- 
luminous per unit gas mass across the starburst region. Column densities derived from both C^^O 
and 1.4 mm dust continuum emission imply that Xco(Maf 2) is about 2-4 times lower than the 
Galactic value, similar to Xco values found for the centers of other large spirals, including our 
own. The weakness of the isotopologues at large galactocentric radii and in the "off-arm" spray 
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regions of Maffei 2, suggest that in these regions either the isotopologues cease to effectively trace 
molecular gas or that the Galactic conversion factor overestimates the molecular column. The lack 
of applicability of the Galactic Xco to the clouds in the center of Maffei 2 is probably due to the 
effect of bar motions and strong tides on the structure and dynamics of these clouds. 

Millimeter continuum emission reveals three prominent locations of star formation with the 
most intense occurring where the molecular bar intersects the nuclear ring. Lyman continuum rates 
of Nlj/c ~ 3-5 X lO^-*^ s~^ are implied for individual regions. The total rate for the entire nucleus is 
Lyc ~ 2.6 X 10^2 g-i^ sFj^ ^0.26 M© yr-^. 

A P-V diagram of the nucleus of Maffei 2 shows a distinct "parallelogram" indicating molecular 
gas response to a barred potential. The morphological and kinematic data confirms Maffei 2 as 
true double barred galaxy. We suggest a bar model where the nuclear gas distribution and velocity 
is governed by a small nuclear bar of r =~110 pc. An upper limit to the mass inflow rate along 
the nuclear bar is dM/dt ^ 0.7 Mq yr^^, enough to drive the current star formation rate seen at 
the end of the bar arms and populate the nuclear ring with gas. The locations of star formation 
and the dense gas in the central region appear to coincide with the location of the xi — X2 orbit 
crossings of the nuclear bar, consistent with dynamical triggering of the the star formation. 

We are grateful to the faculty and staffs of OVRO and BIMA for their support during the 
observations. We thank the referee, Marshall McCall, for a careful and insightful reading of the 
paper. DSM acknowledges support from NSF AST-0506669, the Laboratory for Astronomical 
Imaging at the University of Illinois (NSF AST-0228953) and NRAO. The National Radio As- 
tronomy Observatory is a facility of the National Science Foundation operated under cooperative 
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Table 1. Maffei 2 Basic Data 



Characteristic 


Value 


Reference 


Revised Hubble Class 


SBb(s) pec 


1 


Dynamical Center 


a(J2000) = 02^^41^54^90 ± 0/15 


3 




5( J2000) = +59°36'14."4 ± 2" 






136.5'',-0.3° 


1 


^Isr 


-30 kms~^ 


3 


Adopted Distance 


3.3 Mpc 


4 


Inclination Angle 


67" 


3 


Position Angle 


206" 


3 


M(HI)^ 


4.0 X 10^ Mq 


3 


M(H2)^ 


8.5 X 10^ Mq 


5 



Corrected for adopted distance. 

Corrected for adopted distance and assumed CO conversion factor. 



References. — (1) [Hurt et al.| ( |1993aD ; (2) [Hurt fc Tn^^ ( |1991[ ); (3) 
Hurt, Turner Ho| ( [19961 ); (4) [Fingerhut et al.| ( [20071 ), See text; (5) [Mason 
WilsonI (120041). 
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Table 2. Observational Data 



iransition 




sys 


AT/ 

^ Vchan 


^band 


Beamsize 


Noise 


JJet. 














Level 


b iux 






(K) 


/I — 1 \ 
(km s ) 


(MHz) 


(";°) 


(mK/mJy/Bm) 


(%) 


OVRO: 
















HCN(l-O)'' 


88.63 


300-410 


13.53 


128 


3.8x3.3;-29° 


120/10 


65 


13CO(1-0)^ 


110.20 


230-430 


2.72 


128 


3.9x3.4;-76" 


77/10 


92 


i3C0(2-l)<^ 


220.40 


500-1000 


2.72 


128 


3.3x2.9;-76° 


75/28 


~50 


0^80(1-0)^ 


109.78 


240-430 


10.92 


128 


2.6x2.2;-84" 


130/7.5 




Ci80(2-l)b 


219.56 


300-1000 


5.46 


128 


1.7xl.5;-62° 


150/32 




3.4 mm'^ 


88.92 


300-410 




1000 


2.5x2.5;0° 


17/0.67 




2.7 mm'^ 


109.5 


230-430 




1000 


3.9x3.4;-76° 


3.9/0.50 




1.4 mm^ 


219.3 


300-1000 




1000 


1.7xl.5;-62° 


25/2.5 




BIMA: 
















12CO(1-0)<= 


115.27 


380-1300 


4.07 


172 


3.2x3.1;-14° 


1400/0.15 


105 



Note. — Dates for the observations are ^^CO(l-O), 1994 October 23-1995 January 2; 
^3C0(2-1), 1993 October 26-1994 January 13; 0^^0(1-0) and 0^^0(2-1), 1998 October 19- 
1999 January 5; HCN(l-O), 1999 January 28-1999 March 29; ^^00(1-0), 2004 March 15. 

'^The percentage of the single-dish flux detected by the interferometers. Single-dish inte- 



13 



CO(l-O), ( 


Wehachew, CasoH Sz 


(Wild et al.||l992|. 



''Phase Center #1: Vl5j?=-30 km s'^ a = 02'*38™08^00, <5 = -h59°23'20."0 (B1950), #2: 
a = 02''38'"08^25, <5 = -h59"23'27."0 (B1950). 

^Phase Center #1: Vl5j?=-28 km s''^ a = 02''38'"07^00, 5 = -h59°23'33."0 (B1950) #2: 
a = 02'^38™09^00, <5 = +59°23'40."0 (B1950). 

^Phase Center #1: VL5ij=-28 km s'^ a = 02''38"'07^50, 5 = -h59''23'08."0 (B1950) #2: 
a = 02''38™08^80, 6 = -h59°23'33."0 (B1950) 

''Phase Center #1: VL5ij=-15 km s'^ a = 02'^41™59M9, 5 = -h59°36'46."8 (J2000), #2: 
a = 02'^41"^55^0, S = +59''36'15."0 (J2000), #3: a = 02'^41'"50^31, 6 = +59°35'43."2 
(J2000). 
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Table 3. Measured Intensities 





i2CO(l-0) 

(K km s-i) 


HCN(l-O) 

{K km s-i) 


13CO(1-0) 

{K km s-i) 


i3CO(2-l) 

(K km fi-i) 


C^8O(l-0) 

{K km s-i) 


Ci80(2-1) 

{K km s-'-) 


A 


520 ± 50 


<8.3 


8.5 ±0.9 




< 5.4 




B 


540 ± 50 


26 ±3 


38 ±4 


~ 4.8 


~ 8.7 


~ 5.2 


C 


780 ± 80 


56 ±6 


56 ±6 


9.3 ± 1.8 


23 ±2 


13 ±3 


D 


1300 ± 100 


110 ± 10 


93 ±9 


73± 10 


26 ±3 


18 ±4 


E 


1200 ± 100 


100 ± 10 


120 ± 10 


83 ±20 


28 ±3 


21 ±4 


F 


680 ± 70 


83 ±8 


85 ±9 


41 ±8 


21 ±2 


17 ±3 


G 


530 ± 50 


23 ±2 


63 ±6 


11 ±2 


17±2 


<3.5 


H 


170 ±20 


< 8.3 


<6.8 




< 5.4 





Note. — Uncertainties are the larger of the map uncertainty or the absolute cal- 
ibration uncertainties (assuming to be 10 % for the 3 mm lines and 20 % for the 1 
mm lines). Refer to Figure |2] or Table |4] for the GMC positions. Values are obtained 
from maps convolved to the ^^CO(l-O) beamsize. For these lower resolution maps the 
GMCs are sampled at the "1" component of each GMC (except G which is sampled 
at G3), which corresponds approximately to the peak in the lower resolution maps. 
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Table 4. Giant Molecular Clouds 



GMC 


a, 6 


Q X 6; pa 






rn 

J- pk 


M ■ 




(02''41™;59°36') 


(vc X OC'°) 




{km s^^) 


(K) 


flO*" Mn) 


Al 


57.17,42.4 


51 X 36;46 


78±3 


-66±1 


13 


34 


A2 


57.86,42.2 


29 X 29;120 


82±8 


-59±3 


9 


30 


A3 


57.51,41.6 


52 X 31;23 


95±7 


-67±3 


21 


47 


Bl 


55.95,26.9 


76 X 42;180 


75±2 


-89±1 


9 


42 


B2 


55.75,29.6 


55 X 41;18 


41±2 


-80±1 


17 


10 


C 


55.48,24.8 


110 X 30;26 


61±1 


-84±1 


18 


29 


Dl 


55.14,20.7 


76 X 33;28 


84±1 


-74±1 


31 


46 


D2 


54.78,24.7 


39 X 23;34 


8.6±3 


-54±1 


9 


0.29 


E 


55.08,18.4 


110 X 26;160'' 


llOibl 


-79±1 


19 


87 


F 


54.86,09.6 


95 X 33;36 


51±1 


-2. Oil 


24 


19 


Gl 


54.39,03.2 


65 X 29;36 


47±1 


20±1 


17 


12 


G2 


54.34,00.1 


44x < 17;55 


69±4 


38±2 


9 


<18 


G3 


54.12,59.6^ 


69 X 29;73 


66±2 


38±2 


9 


25 


G4 


53.88,57.9^^ 


< 19x < 17;- 


79±7 


30±3 


12 


<14 


G5 


53.77,00.7 


73 X 22;35 


64±2 


23±1 


14 


22 


HI 


53.76,04.7 


64 X 23;57 


18±1 


14±1 


12 


1.7 


H2 


53.52,01.4 


63 X <17;150 


70±3 


24±2 


9 


<21 



Note. — Positions are based on fitting the uniformly weighted CO(l-O) data. Refer 
to Figure |2] for the locations of each GMC. Uncertainties are la from the least-squares 
Gaussian fits to the data. A GMC is considered unresolved if its deconvolved size is less 
than 1/2 of the beam minor axis. 

aQ2/i4im. 59035' 

'^Size is uncertain due to blending. 
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Table 5. Millimeter/Radio Continuum Flux and Spectral Indices 



Parameter I II III IV V Total 



(ao,5o) 


55.09;21.0 


55.10;17.6 


55.14;16.0 


55.17;13.5 


54.54;01.8 




S-^jj, (4.885) 


14.3±0.7 


19.4±1 


15.2±0.8 


11.4±0.6 


2.7±0.1 


97±5 


SmJy{lA.M) 


7.0±0.4 


9.6±0.5 


7.0±0.4 


5.6±0.3 


1.0±0.2 


48±2 


5„jj,(89.96)'^ 


5.8±0.9 


5.6±0.9 


^3.9 


5.0±0.8 


1.5±0.6 


32±5 


5™j^(110.2)b 


7.0±1 


~7.0^ 


4.7±0.7 


4.1±0.7 


2.7±0.5 


26±3 


5^j2,(219.3) 


19±4 


21±4 


10±3 


~5.3 


9.7±3 


79±16 


al 

«3.4 
"1.4 


-0.6±0.1 
-O.liO.l 
1.5±0.5 


-0.6±0.1 
-0.3±0.1 
~ 1.6 


-0.7±0.1 
< -0.33 
1.1±0.6 


-0.6±0.1 
-O.liO.l 
~ 0.4 


-2.0±0.2 
0.2±0.3 
1.9±0.6 


-O.GzbO.l 
-0.2±0.1 
2.0±0.5 



Note. — Uncertainties are based on the larger of the map noise or 5% for 6 and 2 
cm (Turner & Ho 1994), 15 % absolute calibration at 3 mm (larger uncertainties reflect 
possible errors in line removal), and 20 % at 1 mm. Upper limits are 2a. 

'^Continuum emission is corrected for HCN(l-O) k, HCO"'"(1-0) line flux. 

'^Continuum emission is corrected for ^^CO(l-O) line flux. 

^Uncertain due to confusion with I. 
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Table 6. Star Formation Rates 



Source 




NLyc 


SFR 


^9 7 

^mjy 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


I 


4±0.6 


5±0.6 


0.05±0.01 


15±3 


4.8±0.8 


II 


4±0.6 


5±0.6 


0.05±0.01 


19±4 


e.iibi.o 


III 


2±0.6 


3±0.6 


0.03±0.01 


7±3 


2.3±0.8 


IV 


3±0.6 


4±0.6 


0.04±0.01 


< 2.5 


< 0.81 


V 


<1.2 


<1.6 


<0.02 


10±3 


3.3±0.8 


Total 


20±3 


26±3 


0.26±0.05 


50±10 


16±2.6 



Note. — Table columns: (1) The radio continuum sources 
listed in Table g (2) Free-free flux, ffSmjy, at 88.92 GHz, 
is derived from the spectral flts (Figure [4]). (3) N^^c is the 
number of ionizing photons in units of 10^^ s~^, or ~ 100 07 
stars ( Vacca, Garmany k, Shull|l996 | implied by Smjy (4) 
SFR is the star formation rate in units of Mq yr~^, based 
on the conversion from f^iyc to SFR of Kennicutt (1998). (5) 
Thermal dust flux, ^Smjy at 219.3 GHz is derived from the 
spectral flts (Figure [4]). (6) ^Mjj^ is the total mass in units 
of 10^ Mq implied from the dust flux, ^Smjy, assuming a 
gas-to-dust ratio of 100. 
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Table 7. Observed Line Ratios 



GMC 


12CO(1-0) 
i3CO(l-0) 


12CO(1-0) 
CisO(l-O) 


13CO(1-0) 
Ci«O(l-0) 


12CO(1-0) 
HCN(l-O) 


13C0(2-1) 
lacO(l-O) 


Ci80(2-1) 
CisO(l-O) 


A 


62±20 


> 96 


> 1.6 


> 63 






B 


14±2 


~ 95 


~ 6.7 


21±4 


~0.24 


~0.90 


C 


14±2 


34±7 


2.5±0.4 


14±2 


0.18±0.05 


0.56±0.1 


D 


14±2 


50±6 


3.7±0.5 


12ib2 


0.77±0.2 


0.68±0.2 


E 


lOzbl 


41±4 


4.0±0.5 


12±2 


0.75±0.2 


0.78±0.2 


F 


7.8±1 


32±5 


4.1±0.7 


S.libl 


0.51±0.1 


0.79±0.2 


G 


8.4±1 


31±6 


3.8±0.9 


23±5 


0.28±0.07 


<0.82 


H 


~ 25 


> 31 


> 1.3 


> 20 







Note. — Based on the resolution of the ^'^00(1-0) data. No corrections for differ- 
ences in resolved out flux have been included. 
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Table 9. Maffei 2 Bar Model Parameters 



Characteristic 


Value 


position angle 


29° 


i 


67° 


a a 
^main 


-17° 


^main 


40 km s~^ kpc~^ 


^nuc 


135 km s^^ kpc~^ 


ft b 


7° 


^nuc 




main„ c 
' max 


3.6 kpc 


main„, c 


172 km s"^ 


nuc c 
' max 


155 pc 


nuc c 
Umax 


90 km s~i 


'-'■mam 


1.0 


n ^ 


3.25 


f ■ 

^matn 


0.125 


f 

^nuc 


0.075 




0.125 




0.1875 



^Angle between the main bar and the major axis. 

Angle between the main bar and the nuclear bar. 

'^Brant model parameters for the two axisymmetric 
potentials. The bar core radii are \/2 smaller. 



^Strength of each bar (see Wada 1994 1. 



'^Magnitude of the radial (A) and azimuthal (//) 
damping term (see Wada|1994{ Maciejewski &; Sparke 
19971). 



-37- 




Fig. 1. — Integrated intensity of CO isotopologues in Maffei 2. a) CO(l-O), with contour levels of 
10.0 Jy/beam km s^^ (70 K km s-^ ~ 3a). h) 13CO(1-0), with contour levels of 0.70 Jy/beam 
km s~^ (5.4 K km s^^, 3(t). c) C^^O(l-O), with contour levels of 0.70 Jy/beam km s^^ (5.4 K 
km ~ Scj). d) HCN(l-O) with contour levels of 0.70 Jy/beam km s'^ (8.3 K km s'^, ~ 2a). 
e) 13C0(2-1) with contour levels of 2.5 Jy/beam km s-^ (4.8 K km s'^, ~ 3a). f) 0^^0(2-1) with 
contour levels of 2.5 Jy/beam km s^^ (4.9 K km s^^, 2a). All transitions have been convolved 
to the resolution of the ^^CO(l-O) transitions (3. "9 x 3. "4), shown in panel a). 
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Fig. 2. — a) Uniformly weighted integrated intensity maps of CO(l-O) in Maffei 2. The grayscale 
is in square root stretch ranging from 4.0 Jy bm^^ km s^^ (95 K km s^^) to 60 Jy bm^^ km 
s"^ (1400 K km s"^), for a 2."l x l."9 beam, b) The CO(l-O) map overlaid on the HST F814W 
image of the nucleus. Only contours 118 K km s^^ x (1, 4, 8, 12, 16) are shown. The cross marks 
the location of the dynamical center (Table [T]). c) The CO(l-O) map with locations of the fitted 
GMCs (Table [4]) labeled. Contours are the same as a), d) A zoom in on the central region of the 
molecular bar. CO(l-O) is in contours (steps of 118 K km s~^) overlaid on the 89 GHz continuum 
image (grayscale). Radio continuum (2 cm) source identifications of Tsai et al. (20061 are labeled 
with the symbols given in the legend. Symbol sizes are proportional to the 6 cm peak flux (1 mJy 
bm~^ = 1"; Tsai et al. 



2006). The green cross marks the location of the dynamical center. 




Fig. 3. — Radio and millimeter continuum maps of the nucleus of Maffei 2. Contours are ±2"/^, 
n=0,l, 2... times the given lowest contour level. Dashed lines are negative contours, a) 2 cm 
map of Turner &; Ho (1994), with the lowest contour at 0.45 mJy beam^^ (2<7)- Roman numerals 
identify the five major continuum sources, b) The 89 GHz (3.4 mm) map contoured in steps of 
1.2 mJy beam~^ (2o")- Line contamination from HCN(l-O) and HCO^(l-O) has been removed, c) 
The 110 GHz (2.7 mm) map contoured in steps of 1.0 mJy beam~^ (2o"). Line contamination from 
^■^00(1-0) has been removed, d) The 219 GHz (1.4 mm) map convolved to the resolution of the 
110 GHz map. The lowest contour is 6.0 mJy beam~^ (2<7)- Beamsizes are indicated at the bottom 
left of each map, and in Table |2} 
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Fig. 4. — Spectral energy distributions for the millimeter continuum radio sources in Maffei 2. 
Data are fit with synchrotron (dashed line, a = —0.7), free-free (dotted line, a = —0.1) and dust 
(dot-dashed line, a = -1-3.5) emission components (see text). The cross in Source III represents an 
upper limit. Triangles mark the contribution of the compact emission towards each source from 
higher resolution images (Tsai et al. 2006). 
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Fig. 5. — CO and HCN line ratios in Maffei 2. For comparison the -'^■^CO(l-O) integrated intensity 
is overlaid on all planes in linear contours of 2.0 Jy/beam km s~^ (15 K km s~^). The resolution 
of all the ratio maps is 3. "9 x 3." 4. The greyscale range of the ratio is noted in the wedge at the 
top of each panel, (a) CO(1-0)/13CO(1-0). (b) CO(l-0)/Ci8O(l-0). (c) CO(1-0)/HCN(1-0). (d) 
1=^CO(2-1)/13CO(1-0). (e) Ci«O(2-l)/Ci«O(l-0). (f) i3CO(l-0)/Ci«O(l-0). 
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Fig. 6. — Average peak T^b ratios for selected regions along the central ring of the P-V diagram. 
The schematic representation of the P-V diagram has the same geometry as Figure 10 The letters 
in parentheses in the upper left label the primary molecular peak in each region. 
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02 41 58 56 54 52 02 41 58 56 54 52^ 02 41 58 56 54 52 



Fig. 7. — The blueshifted half of Maffei 2's ^^CO(l-O) channel maps. Contours are multiples of 30 
mJy/beam ( 2a). The beamsize is given in the bottom left of the first plane. 
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Fig. 8. — The redshifted half of Maffei 2's ^^CO(l-O) channel maps. Contours given in Figure 
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Fig. 9. — Maffei 2 HCN(l-O) channel maps. Contours are in increments of the 2a value, 20 
mJy/beam. The beamsize is given in the bottom left of the first plane. 
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Fig. 10.— The a) CO(l-O), h) 13CO(1-0) and c) HCN(l-O) Position- Velocity (P-V) diagrams taken 
along the major axis of Maffei 2. The zero velocity corresponds to -20.5 km s~^ (LSR). The zero 
position corresponds to 02:41:55; 59:36:10. Northeast along the galaxy is at the top of the figure. 
Contours are ~ 2cr. Labels in identify the location of each GMC in position and velocity. 
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Fig. 11. — An epicyclic double bar model for Maffei 2. See ^6.1 for a description and Table [9] 
for parameters, a) The 2MASS infrared K band image of Maffei 2. h) Nuclear molecular gas 
morphology traced by CO(l-O), showing only the 118 (xl,4,8, 12) K km s~^ contours, overlaid off 
the 814W HST image, cj A model of the inner portion of the large scale bar on the same scale as 
a). Predicted positions of the gas peaks follow the locations of highest point density. Note that 
the northeastern NIR arm is tidally disturbed, and as a result the model agreement is poorer, d) 
A zoom in on the nuclear bar portion of the model, As in c) except displayed face-on. f) As in 
d) except viewed face-on. 
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Fig. 12. — The nuclear velocity field of Maffei 2. a) The ^^CO(l-O) position-velocity diagram taken 
along the major axis (Figure [Toja) . h) The predicted P-V diagram for the nuclear bar on the same 
scale as a), showing the well known 'parallelogram' characteristic of barred motion. Color coding 
shows the portion of the velocity field associated with each of the components labeled in cj. c) 
Regions of high point density correspond to predicted (and observed) regions of gas concentration. 
The three main regions of predicted gas concentration are color coded — the northern and southern 
arms are green and blue, respectively and the central ring is pink. 
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Fig. 13. — Rotation curve of MafFei 2. a) The model rotation curve (thick red line) and observed 
HI rotation curve (red errorbars; Hurt, Turner &; Ho 1996) is displayed. Blue (yellow) dot-dashed 
lines correspond to the pattern speeds of the nuclear (main) bar. Associated resonance curves 
are labeled, b) A plot of dynamical properties as a function of galactocentric radius based on the 
modeled rotation curve. Shown are the azimuthally averaged observed velocity dispersion (green), 
molecular gas surface density (dark blue) from ^^CO(l-O) assuming Tg^ — 5 K and an t — 1 (see 
text), epicyclic frequency (pink), implied enclosed dynamical mass (light blue) and Toomre's Q (see 
text) (red). Units for each curve is given in the figure. 




Fig. 14. — LVG model solutions for six locations across the nucleus. Abundances per velocity 
gradient, Xco/dv/dr^ between 10^^ '^ — 10^^ have been modeled. ^Xco — 8 x 10^^ for [CO] / [H2], 
(1/60) (8 X 10-5) for [i^CO] / [H2], (1/250) (8 x 10"^) for [C^^O] / [H2], and a velocity gradient of 
1 km s~^ pc"^ are assumed (^7.1.1). Solutions are based on the ^^CO(l-O) resolution. Line ratios 
have been corrected for resolved-out flux. Contours represent the la confidence solutions. Red 
contours display the isotopic line ratios (thick: CO(l-0)/i3co(l-0); thin: CO(l-0)/Ci8O(l-0)) 
and blue contours the AJ line ratios (thick: 13CO(2-1)/13CO(1-0); thin: Ci8O(2-l)/Ci*^O(l-0)). 
Green contours display the allowed parameter space based on the ^'^CO(1-0)/HCN(1-0) line ratio, 
assuming the filling factor of the HCN(l-O) emission is equal to that of the ^^CO and C^^^O emission. 
Xhcn / dv / dr = 2.1 x 10~®, is assumed, consistent with a Galactic Center HCN abundance and a 
velocity gradient of ~ 1 km s~^ pc~^ (eg., Paglione et al.|l998 |. The yellow line marks the CO(l-O) 
peak brightness temperature assuming a filling factor of unity for the uniformly weighted beamsize. 
Crosses mark the adopted best fit solutions. 



